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We compute the spectral density in the normal phase of an interacting homogenous Fermi gas 
using a T-matrix approximation. We fit the quasiparticle peaks of the spectral density to BCS-like 
dispersion relations, and extract estimates of a “pseudo-gap” energy scale and an effective Fermi- 
wavevector as a function of interaction strength. We find that the effective Fermi-wavevector of 
the quasiparticles vanishes when the inverse scattering length exceeds some positive threshold. We 
also find that near unitarity the quasiparticle lifetimes, estimated from the widths of the peaks in 
the spectral density, approach values on the order of the inverse Fermi-energy. These results are 
consistent with the “breakdown of Fermi liquid theory” observed in recent experiments. 
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I. INTRODUCTION 

Interacting degenerate Fermi gases have attracted con¬ 
tinued interest since their experimental realization over 
a decade ago [l]. Below a critical temperature T c these 
gases exhibit a superfluid state which can be continuously 
tuned via Feshbach resonances [|!] from a BCS state of 
Cooper pairs, to a BEC state of tightly bound bosonic 
molecules @]. While the physics of the superfluid phase 
in this BCS-BEC crossover is well established Q, much 
less is known about the normal phase at temperatures 
above T c . 

One theme that has emerged from theoretical and ex¬ 
perimental investigations at T > T c is the idea of a “pseu¬ 
dogap” phase in the middle of the crossover, where the 
density of states at the Fermi energy is suppressed due 
to strong many-body pairing effects. If and how this 
phase emerges in the BEC-BCS crossover has long been a 
source of experimental and theoretical investigation (see 
Ref. |5j for a comprehensive review). More recently, there 
has been interest pfToi in a related question: How does 
Fermi-liquid theory [llj, which is expected to be valid in 
in the BCS regime, break down when crossing over to the 
BEC regime where the normal phase is a gas of weakly in¬ 
teracting bosons? There has been some disagreement in 
the conclusions drawn from experiments which use radio¬ 
frequency (RF) spectroscopy to probe the single particle 
spectral density in harmonically trapped systems. The 
observed RF spectrum in the crossover region seems to 
be well described by both Fermi-liquid theory Q and 
theories displaying a pseudogap phase fb, @]- A more 
recent experiment [9] probed the spectral density in a 
nearly homogenous system which avoids the density in¬ 
homogeneity of the trapped systems that can obscure 
features in the RF spectrum [121 • This experiment found 
evidence that the well defined quasi-particles one expects 
from Fermi-liquid theory become absent as the interac¬ 
tions are tuned from the BCS to the BEC side. 

In this paper we address this problem theoretically 
by computing the single particle spectral density within 
a T-matrix approximation as a function of interaction 


strength. We find near unitarity that the quasiparticles 
represented by peaks in the spectral density have short 
lifetimes at the Fermi-wavevector. We also find that an 
effective Fermi-wavevector extracted from the shape of 
the quasiparticle dispersions vanishes when one moves 
sufficiently deep into the BEC regime. Both of these ob¬ 
servations point toward a breakdown in the Fermi-liquid 
description of the normal phase in which there is a well- 
defined Fermi surface and long-lived quasiparticles at the 
Fermi-wavevector. Previous theoretical works using re¬ 
lated T-matrix approximations have similarly discussed 
the shape and widths of peaks in the spectral density 1131- 
[Tl . Much of that work focused on the BEC regime [ll| 
or on the temperature dependence of the spectral den¬ 
sity at a few discrete values of the interaction strength 
[14 ilT . We extend these results and systematically ex¬ 
plore the dependence of the spectral density on interac¬ 
tion strength. 

This paper is organized as follows. In Sec. [IT] we discuss 
the T-matrix approximation in detail. In Sec. lIIIl we show 
numerical results for the spectral densities and describe 
our procedure for analyzing the quasiparticle dispersions 
and lifetimes. In Sec. m we discuss our results and con¬ 
clude. 


II. T-MATRIX APPROXIMATION 

In this paper we consider a Hamiltonian that models 
two component fermions with attractive s-wave interac¬ 
tions 

H = X e kaL°k<T + X a kt a pt a P- t Fl' Ct k+qt (1) 

k kpq 
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where £k = ^7 — U and the interaction strength g is re¬ 
lated to the scattering length a by g -1 = — -I pf- 

(we set H = 1 and Boltzmann’s constant ks = 1 through¬ 
out). 

We use a non-self-consistent T-matrix approximation 
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FIG. 1: Diagrammatic representation of the T-matrix approximation used in this paper (Eq. (|2|)). 


, to calculate the single-particle self energy E: 


E(q, iuj n ) = ^ E E k ( k ’ *Hrn)Go(q - k, i^m - iu n ) 
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where u) n = (2n+l)7rT and Q m = 2rmrT (with integers m 
and n) are fermionic and bosonic Matsubara frequencies 
at temperature T. Go(k, iui n ) = l/(iui n — e k ) denotes the 
bare single particle propagator. The density is given by 

n ~ = 2 v EE G ( k + w ") ( 3 ) 

k n 


where G 1 = G 0 1 — E. The critical temperature is set 
by a pairing instability condition (the Thouless criterion) 


T-matrix theories [l§| (where the bare propagators in 
Eq. © are replaced by the fully dressed propagator G) 
produce results with similar qualitative features. Since 
we are mainly concerned with discussing the qualitative 
physics within this regime, we use the relatively simple 
approximation described by Eq. ©. 

After performing the Matsubara sums in Eq. © and 
taking iu n —> lining ui + iS we arrive at the following 
expressions for the imaginary part of the self energy: 

ImE(q,w) = (5) 

/ d 3 k 

( 2 ? r) 3 ^( k + <+ w + £ k )[/b(w + e k ) + /f(c k)] 

where f F (e) = (exp(e/T) + 1) _1 , / B (e) = (exp(e/T) - 
1) _1 , and 

=ImT(q,a;) = (6) 
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The expressions in Eq. © are derived by summing the 
infinite subset of Feynman diagrams that include scatter¬ 
ing processes occurring in the vacuum two-body problem 
(this shown in Fig. ©. This approximation, which ne¬ 
glects interactions between pairs of fermions, neverthe¬ 
less accurately models the physics in the weakly interact¬ 
ing BCS regime << — 1, where the normal phase 
is a Fermi-liquid) and in the strongly interacting BEC 
regime (^+^ >> 1, where the normal phase near T c is a 
gas of free bosonic molecules of mass 2m). In particular 
in the deep BCS regime, the Thouless criterion is equiv¬ 
alent to the expression for T c from BCS theory; in the 
deep BEC regime, the Thouless criterion together with 
Eq. © yield the correct condensation temperature for a 
gas of non-interacting bosons [la¬ 
in the crossover regime (—1 < < 1) this approx¬ 

imation is less well-controlled and there is no a priori 
reason to expect it to be accurate. However previous 
studies using ab initio techniques [l8[ and self-consistent 


with d —>■ 0. We numerically evaluate the integrals in 
Eqs. © and © and extract the real part of the retarded 
self energy using a Kramers-Kronig relation. The single 
particle spectral density A(k, ui) is then given by 

A(k, ui) = — 2ImG(k, to) 

—2ImE(k, u>) (7) 

\lu — e k — ReE(k, w)] 2 + [ImE(k, w)] 2 

III. RESULTS 
A. Spectral Density 

Fig. [2] shows the spectral density A(k, E) calculated 
for interaction strengths —0.5 < —< 0.5 at temper¬ 
atures T = 1.1T C . This choice of temperature was mo¬ 
tivated by recent experiments performed slightly above 
T c 0. In the weakly interacting BCS limit (for instance, 
the bottom right panel in Fig, © we observe a strong 











3 



FIG. 2: (Color online) Spectral density in the normal phase for various interaction strengths at temperature T = 1.1T C . Lighter 
colors correspond to higher spectral density. 


quadratically dispersing quasiparticle peak along with a 
faint peak at E > 0, k < kF- 

As the interactions become stronger we see two clearly 
identifiable quasiparticle dispersions, one with E > 0 and 
another with E < 0. These two branches are sep arated 
in energy by a depression (or “pseudogap” [H], |2IJ |2l| ) in 
the density of states centered at E = 0, k = Uf- Finally, 
for very strong interactions (the top panels in Fig. [2]) the 
two quasiparticle branches become further separated in 
energy and the spectral density is nearly fully suppressed 
near E = 0, k = Uf- Similar results calculated from 
various T-matrix approximations are shown in Refs. 

MM- 

In the BCS regime there is a simple cartoon picture [22j 
of these results: One can add a fermion by either occu¬ 
pying a “normal” fermion particle state with dispersion 
E « ek or by creating a “pair” of energy 2? pa ir(p) and 
annihilating a fermion with energy e p _k- The latter exci¬ 
tations are broad as the pairs will be created with a range 
of p. In this regime, the pair has vanishingly small energy 
-Epair(O) ~ 0, so we expect the “pair” peak to roughly 
track E ss — e k . As the interactions are increased, there is 
hybridization between these two branches that produces 
a psuedogap analogous to the gap that is produced by 
the hybridization of particle and hole states in the super- 
fluid phase. Finally in the BEC regime, E pa i r < 0, which 
shifts the branch associated with pairs entirely to nega¬ 
tive energies; the positive energy branch in this regime 
corresponds to adding one “normal” fermion to a sea of 
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FIG. 3: (Color online) Peak locations of the spectral density 
at 1/kFd = 0.18 (T = l.lTo). The solid lines show least- 
squares fits to the dispersion relation in Eq. ©. The board 
and asymmetric nature of the peaks prevents us from deter¬ 
mining the lower branch dispersion in the range 1.0 kF < k < 
1.25 k F . 


paired fermions for all k. 
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B. Quasiparticle Dispersions and Lifetimes 

To describe the qualitative features of the quasiparticle 
dispersions, we fit the peaks in the spectral density to 
a phenomenological model obtained from the dispersion 
relations expected in the superfluid phase: 


£peak(*0± = (8) 

where E pea k(k)± are the peak locations as a func¬ 
tion of momentum for the positive and negative energy 
branches. A± parametrizes the magnitude of the psuedo- 
gap energy scale, and kL ± parametrizes the location of 
the effective Fermi-wavevector. Note that Ul (called the 
“Luttinger-wavevector” in Refs. @, 0 ) may not equal 
k F as defined in Eq. ®. An example of such a fit is 
shown in Fig. [3] (where 1 /kpa = 0.18, A_ = 0.95-Ef, 
A + = 0.94 Ep, k F _ = 1.01 kp, and = 0.34 A:f). 

Fits like these have been performed in previous work 
in the BEC regime [l3j and for a few select values of 
the interaction strength at different temperatures |l5j . 
Here we focus on the dependence of the fit parameters 
on interaction strength and perform fits for a number of 
different values of 1/fcpa in the crossover regime. Simi¬ 
lar strategies have been used to capture the qualitative 
features of experimental data for harmonically trapped 
fermions [23| • 

Figs. S] and [5] show the fitting parameters A± and k F± -, 
respectively, as a function of interaction strength. Each 
point was exacted by performing a least squares mini¬ 
mization of the difference between Eq. ® and the loca¬ 
tion of the quasiparticle peaks in the spectral densities 
in Fig. [2] The psuedogap scale is maximal on the BEC 
side of resonance, smoothly dropping as one approaches 
the BCS side. These parameters are determined by a 
global fit, so they are sensitive not only to features near 
the Fermi wave-vector, but to all k. While k F+ and 
are generically quite different, we find A + ps A_. 

To estimate the lifetimes of the quasiparticles we ex¬ 
tract the widths of the spectral peaks at k = kp by fitting 
the peaks in the function A{k = kp,E) to Lorentzian 
functions. Fig. [3 shows an example of such a fit for 
—= 0.18. Fig. [6] shows the quasiparticle widths r” 1 
as a function of interaction strength. We note that for 
< 0 the negative energy branch cannot clearly be 
separated from the positive energy branch at k = for 
these interaction strengths we only plot one value of the 
quasiparticle lifetime. 


IV. DISCUSSION AND CONCLUSIONS 
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FIG. 4: (Color online) Pseudogap energy scale A_ (orange) 
and A+ (blue) extracted from fits of Eq. ((HJl to the peak lo¬ 
cations in the spectral density (T = 1.1T C ). 
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FIG. 5: (Color online) Effective Fermi-wave vectors (or¬ 
ange) and kL + (blue) extracted from fits of Eq. (|8]) to the 
peak locations in the spectral density (T = 1.1T C ). 
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Our results demonstrate a breakdown of Fermi-liquid 
theory in two key respects. First, as shown in Fig. [6j the 
inverse lifetimes of the quasi-particles at kp increase as 
one approaches unitarity from the BCS-side and even¬ 
tually reaches a maximum of value on the order of Ep. 


FIG. 6: (Color online) Inverse quasi particle lifetimes r _1 as 
a function of interaction strength (T = 1.1T C ). The orange 
(blue) points are the spectral widths of the negative (positive) 
energy branch of the quasiparticle spectral peak at k = k F . 
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FIG. 7: (Color online) A(k = k F ,E) at ^ = 0.18 (T = 
1.1T C ). The solid red line is a fit of the data to the sum of 
two Lorenzians. The two widths extracted from these fits 
provide estimates of the inverse quasi-particle lifetimes r _1 . 
t -1 at other interaction strengths are plotted in Fig. [U] 

Fermi-liquid theory is predicated on a vanishingly small 
inverse lifetime at the Fermi-surface. This breakdown is 
due to both the high temperature and the strong inter¬ 
actions. 

Second, the effective Fermi-wavevector k F , which helps 
identify the presence of a remnant Fermi surface |l5[ , 
vanishes for the upper branch at ~ 0.3 (see Fig. 0. 
We note that the effective Fermi-wavevector for the bot¬ 
tom branch (k F _) remains fixed at k F for the range of 
interaction strengths shown in Fig. [5]). However, other 
studies have shown that at higher interaction strengths 
(iTa ~ O-'O’ vanishes as well [§}. and k Fjr 

are associated with different branches of excitations and 
therefore have distinct physical interpretations. is 
associated with particle excitations, while is associ¬ 
ated with holes. On the BEC side of resonance, one can 
add a fermion without disturbing the pairs, but adding a 
hole requires breaking pairs. Thus many-body effects are 
less important for particles, and fc_L + vanishes in a regime 
where k F _ is still essentially equal to k F . 

A recent experiment |9:j used RF spectroscopy to probe 
the spectral density of a nearly homogenous Fermi gas in 
the normal phase as a function of interaction strength. 
The RF signal reported in this work came primarily from 
the negative energy branch of the single particle excita¬ 
tions since any signal from positive energy excitations is 
suppressed by a Fermi factor. The authors fit their data 
to a two-mode model with a narrow quadratically dis¬ 
persing peak expected from a Fermi-liquid theory (with 


weight Z < 1) and a broad “incoherent background” 
(with weight (1 — Z)) corresponding to the spectral weight 
of weakly interacting bosonic molecules in the normal 
phase. They find that the best fit to their data at in¬ 
teraction strengths > 0.3 has Z « 0. The authors 
conclude that this is a signal of a breakdown in the Fermi- 
liquid description of their data. 

We caution however that because their two mode fit 
was performed on a single broad dispersing peak (com¬ 
ing from the negative energy branch of excitations), it is 
somewhat difficult to interpret the meaning of this van¬ 
ishing Fermi-liquid contribution. Moreover, good fits to 
the data came at the expense of some unrealistic values of 
the fitting parameters. In particular, the parameter cor¬ 
responding to the temperature of the incoherent bosonic 
contribution had best fit values nearly four times greater 
than the estimated temperature of the gas. 

One possible interpretation of this two-mode fitting 
procedure is that Z roughly measures the weight of the 
component of the RF signal which disperses as E ~ k 1 2 3 . 
Assuming that the peaks of the RF data are well de¬ 
scribed by the BCS-like dispersion E pea k(k)- given in 
Eq. (0), the quadratically dispersing component vanishes 
when k F _ —> 0. The vanishing of Z in the two-mode 
model then seems to roughly correspond to the vanishing 
of k L _ that occurs as one approaches the BEC regime. 
We recommend fitting the peaks of the RF data to the 
BCS dispersion relations in Eq. ( 0 ) to directly extract 
kL_ and to compare the experimental results with our 
T-matrix calculations. 

Given the high temperatures and strong interactions 
in the normal state near unitarity, it is not surprising 
that many of the features expected of a Fermi liquid 
are absent. It appears that the T-matrix approximation 
captures this physics. It remains to be seen if an alter¬ 
native framework can replace these Fermi-liquid ideas. 
The most tantalizing steps in that direction come from 
exploring universal bounds on transport coefficients [24j ■ 
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